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Abstract: 



The Q 2 dependence and the total cross sections for charged and neutral current 
processes are measured in e p reactions for transverse momenta of the outgoing 
lepton larger than 25 GeV. Comparable size of cross sections for the neutral 
current process and for the weak charged current process are observed above 
Q 2 « 5000 GeV 2 . Using the shape and magnitude of the charged current cross 
section we determine a propagator mass of mw = 84 GeV. 
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1 Introduction 



The unification of weak and electromagnetic forces is an important ingredient of the Standard 
Model. As a consequence the cross sections for neutral current (NC) and charged current 
(CC) processes 

e ± p — > e ± + X and e ± p +X (1) 

(X denotes the hadronic system recoiling against the lepton) are predicted to be of compa- 
rable size when probing small distances. The HI experiment at the electron-proton collider 
HERA permits the study of these processes at large squared four momentum transfer Q 2 
using the same detector for the measurement of both processes. 

The generalized cross sections for deep-inelastic e ± p-scattering can be written as M: 
d 2 a(e ± p) 2na 2 r, ,9,^ n9n _\ 

Here x and y are the Bjorken scaling variables which are related by Q 2 = sxy. The center 
of mass energy squared s is given by the product of the electron and proton beam energies 
s = 4E e Ep. The electromagnetic coupling constant is denoted a and JF 2 3 are generalized 
structure functions. These cross sections contain not only the proton structure functions 
but also the electroweak coupling constants and propagator terms. The contribution of the 
longitudinal structure function F L is expected to be small at high Q 2 and is neglected here. 

The NC cross section is governed by the photon propagator while the Z° exchange affects 
the cross sections significantly only at Q 2 ~ m| - Also the contribution of the Z° exchange 
enters with different signs for e + p and e~p scattering. The difference between the two NC 
cross sections therefore provides a method for detecting these electroweak effects. However, 
the accumulated statistics of the NC data are not yet sufficient to distinguish the dominant 
photon contribution from Z°— exchange. 

For CC processes the lowest order cross section can be simplified to 



dV£f _ Gl 1 



dxdQ 2 2tt(1 + Q 2 /m 



2 \ 

W) 



dV£f _ G\ 1 



dxdQ 2 2tt (1 + Q 2 /m 



W) 



£ q(x,Q 2 ) + £ (l-y) 2 q(x,Q 2 

9= u > c q=d,s,b 



E ^,Q 2 )+ E (l-y) 2 q(x,Q 2 ) 

q=u,c q=d,s,b 



(3) 



(4) 



where q and q are the quark and anti-quark densities of the proton. G^ is the Fermi constant 
known from the measurement of the muon lifetime. In the propagator term mw denotes the 
mass of the W bosons which mediate the weak force. 

In our previous publications total charged current cross section measurements have been 
presented |2|, f| . In this analysis we study the dependence of the charged current cross section 
on the four momentum transfer. Such an analysis has been presented earlier by the ZEUS 
Collaboration [|J for 1993 e~p. Here we present results obtained from both e~p and e + p 
data recorded in 1993 and 1994. From the cross section and the shape of the Q 2 distribution 
we derive m w . 

Neutral current cross sections at high Q 2 have already been precisely measured using the 
information provided by the scattered lepton only || or together with the hadronic final 
state 0. For the comparison with charged current processes we use here the hadronic final 
state for event selection as well as for kinematic reconstruction. This procedure reduces 
uncertainties arising from different systematic effects. 



2 Experimental conditions 



HERA is an ep-collider where electron or positron beams collide with 820 GeV protons. 
This analysis is based on data taken in 1993 and 1994 using the HI detector. The integrated 
luminosities are 0.33 pb _1 (1993) and 0.36 pb _1 (1994) for e~p running and 2.70 pb _1 (1994) 
for e + p running. The lepton beam energies were 26.7 GeV in 1993 and 27.5 GeV in 1994. 

HI is a multi-purpose detector well optimized to measure large Q 2 events in deep inelastic 
scattering. A detailed description of the HI detector and its performance can be found in 
0. We discuss those aspects which are relevant to the present analysis. 

A system of central and forward drift chambers is used to determine the interaction vertex 
of an event. The chambers measure tracks from charged particles in the polar angular range 
of 7° — 165° with respect to the proton beam direction (which is the ©^-direction). The drift 
chambers are complemented by several layers of proportional chambers which provide fast 
trigger signals. 

The energy and angles of the hadronic final state particles as well as the scattered electron 
are measured in the highly segmented (~ 45000 cells) liquid argon (LAr) calorimeter ||. It 
covers a polar angular range between 4° and 153°. The total depth of the calorimeter 
varies between 4.5 and 8 interaction lengths. The resolution for single pions and electrons 



is a(E)/E » 0.50/y/(E/Ge V) © 0.02 and a(E)/E « 0.11/^(£/GeV) © 0.01 respectively, 
as measured in test beam M. The hadron energy scale is verified within 4% with data 



taken at HERA ||10|| . For trigger purposes the LAr calorimeter is read out in a coarser 
granularity via flash ADCs. They provide a fast measurement of energies and allow triggering 
on missing transverse momenta (CC trigger) or on large electromagnetic energy deposition 
(NC trigger). 

Both the LAr calorimeter and the central drift chambers are surrounded by a super- 
conducting solenoid. The iron yoke which returns the magnetic flux is instrumented with 
streamer tubes and is used as a muon detector. 

A luminosity system measures the rate of the small angle Bremsstrahlung process 
ep — > ep7 |TT| and allows the determination of the luminosity with a precision of 1.5% 
(3.5% in 1993). 



3 Data Samples 

The selection starts from all triggered events and is based on information of the hadronic 
final state both for CC and NC interactions. The presence or absence of a charged scattered 
lepton is used only to classify an event as a candidate for a CC or a NC process. 

Scattered electrons are identified in the LAr calorimeter by the lateral extension and 
energy density of the electromagnetic shower and the requirement that a reconstructed track 
points to the shower. The minimal transverse energy deposition of the electron candidate 



is required to be 5 GeV [12j. The search is restricted to scattering angles of the charged 
lepton d e < 150°, which is well inside the angular acceptance of the liquid argon calorimeter. 
The efficiency of the electron identification is determined by comparing independent electron 
identification algorithms, either based on calorimeter information only or on a track criterion. 



A total electron identification efficiency of 98.7 ± 0.3% is found for the finally analysed NC 
sample. 

After the neutral current candidates are identified all information about the electron in 
the track detectors and calorimeters is ignored in the subsequent analysis. The additional 
selection described below is applied to both samples on the hadronic final state only. The 
selection steps are as follows: 

1. The hadronic transverse momentum sum p^h measures the missing transverse momen- 
tum with respect to the incoming particle direction in each event and is used as a main 
selection cut requiring 



PT,k = /(E^) 2 + (E^) 2 > 25 GeV . (5) 

y i i 

The computation uses the event vertex position and energies and positions of LAr 
calorimeter cells (index i) of the hadronic final state. Neglecting particle masses the 
three components of the momentum vector are then defined by {p x ,Py,Pz)- 

2. For both event types the CC trigger condition has to be fulfilled. For CC events this 
signal is derived from the LAr calorimeter trigger energies, for NC events it is evaluated 
from the same measurement after electron subtraction. To assign the events to a 
specific crossing of electron and proton bunches, a trigger signal from the proportional 
chambers is required. 

3. A vertex is reconstructed from measured tracks. This vertex position has to be within 
a distance of 35 cm from the nominal interaction point along the beam axis. 

4. High energy cosmic or proton beam related muons may cause large energy deposits in 
the calorimeter and thus fulfill the selection criteria. Based on the topology of these 
deposits and event timing with respect to the beam crossing background filters have 
been developed which remove these events efficiently [[3]. 

The efficiency for these four selection criteria is determined from the initial NC sample. The 
selection efficiencies are found to be between 80% and 85% for the four processes with an 
uncertainty of 3%. All selection efficiencies depend only little on kinematic properties of the 
events. The final event statistics is given in Tab. [L| 

The CC-candidate sample is subject to a visual inspection by which residual background 
from incoming muons and from misclassified NC events is identified and removed [Q. Pho- 
toproduction events 7p — > hadrons in which the scattered electron remains undetected 
(Q 2 ~ 0) are a potential background for the CC sample if the hadronic final state is badly 
reconstructed due to particle losses in the beam pipe or resolution effects. This background 
has been estimated from data and was found to be less than 0.5% as described in [||, [12]. 
Therefore it has been neglected in further analysis here. 

The total background rate for NC events is 2.0 ± 0.6% and has been statistically sub- 
tracted from the sample. The largest background contribution is due to 7p events with many 
hadrons in the final state, one of which is misidentified as an electron. For some of these 
events the scattered electron is detected in the luminosity system of the HI detector. From a 
study of these tagged events we find a background fraction of 1.4 ± 0.6% []P2| . Events of the 



type e ± p — > e e + e X have been found by a visual inspection of the sample. They contribute 
an additional 0.3 ±0.2% to the background. The background from incoming muons amounts 
to 0.3 ±0.1%. 



4 Data Analysis 



The final state kinematics is fully specified by two independent variables. Following the 
method proposed in [Tj| we use the transverse momentum sum pj-,h (cf. eq. (^)) and the 
energy transfer y^. 

y h = J2(Ei-p Zi )/2E e . (6) 



Here also particle masses are neglected and the summation is performed over all LAr 
calorimeter cells i of the hadronic final state. 

For comparison to data and for studies of migration effects and unfolding we use a detailed 
detector simulation. The electroweak processes are generated by the DJANGO Monte Carlo 
program that includes QED radiative effects of order 0(a) jnj] with proton parton densities 
of the MRSH parameterization [15|. The program makes use of the color dipole model as 



implemented in ARIADNE [IB| for hadronization. 

In Fig. 1 the measurements of pr,h and yh are compared with the expectations from the 
simulation for the NC sample. The simulation agrees well with the measured distributions 
within the statistical errors. 
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Figure 1: Neutral current e + p data (filled circles) are compared to a Monte Carlo (MC) simulation 
(histogram) normalized to the observed number of events. Shown is the transverse momentum sum 
PT,h of the hadronic final state is shown in (a) and yh in (b) . For the data only statistical errors 
are given. The statistical errors from the MC simulation are negligible. 

The squared four momentum transfer Q 2 is related to p\ h and y^ by 

Ql = J^- (7) 

l -yh 

The relative resolution for Ql is about 20%. For experimental cross checks the measurement 
of Q\ by the scattered electron in NC events can be compared with Q\ determined from the 



hadronic final state. The lepton scattering angle i? e and energy E' e are used for calculating 
the kinematics with: Q 2 = 2E e E' e (l + cos($ e )). The Q 2 e measurement has a resolution 
AQl/Ql pa 3% . In Fig. 2 we compare the measured ratio Q 2 h /Q 2 with the Monte Carlo 
simulation. The comparison shows that the hadronic measurement of the kinematics is well 
described by the simulation. The distribution is not centered around one due to losses in 
transverse momentum of the hadronic final state into the beam pipe. 
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Figure 2: The ratio of the squared four momentum transfer as determined from the hadronic final 
state Q\ and from the scattered electron Q 2 for NC e + p events compared to the Monte Carlo 
prediction. For the data only statistical errors are given. The statistical errors from the MC 
simulation are negligible. 



We apply a matrix unfolding technique [17] to determine the true Q 2 distribution from the 
measured Q 2 m distribution. The true squared four momentum Q 2 is defined by the outgoing 
lepton. For the unfolding we make use of the relation 



dQ 2 A(Q 2 m ,Q 2 



da(Q 2 
] dQ 2 



^(Qmy Q 2 ) i s t ne detector response function averaged over the variable x. Consequently it 
has to be determined separately for the different processes under study. For finite Q 2 bins 
equation (pi) can be converted into a sum: 



(9) 



with bin indices fi and v and the transfer matrix T^ v . The transfer matrix elements are 
determined by Monte Carlo simulation. They contain the transformation of the cut in 
PT,h into the corresponding experiment independent cut in the true pt of the scattered 
lepton. The calculated transfer matrix is cross checked for NC data using the complementary 
information provided by the electron. The binning was chosen such (cf. Tab. Ef), that the 



majority (typically 60%) of the events are reconstructed into the original Q 2 bin. Inverting 
relation @ and applying it to the measured binned distribution gives the true distribution 
together with its covariance matrix. 

In binned distributions the information on individual events is lost. To exploit the full 
information on my/ contained in the shape of the CC Q 2 distribution and the total cross 



section we use the method of extended maximum likelihood UI8| . The likelihood function is 
defined as: 

- P (0 2 ) e~ M ■ Af n 

L = TT^'- f- (10) 



where i runs over n individual events. For an event i the distribution Pj(Q^) in four mo- 
mentum Q 2 m is given by the convolution of the differential cross section da(Q 2 ,mw)/dQ 2 
and the detector response function AiQ^Q 2 ) (cf. eq. (||)). The predicted number of events 
for a propagator mass mw and an integrated luminosity £ is Af = £ J P(Q 2 n )dQ 2 n . The 
propagator mass is then extracted from the data by maximizing the likelihood function. 



5 Results 

In Fig. 3 and Tab. [I] we present the result of the unfolding of the Q 2 -differential cross 
sections for the e^p data for NC and CC for px > 25 GeV. The theoretical prediction 
was obtained with the HERACLES |1J generator using the MRSA' |TJj proton parton 



densities. The uncertainty of the prediction is about 4% due to the imperfect knowledge 
of the parton densities. Note that the data are not corrected for radiative effects. Born 
cross sections can be obtained using daBom/dQ 2 = da/dQ 2 {l + <5b ™) (cf. Tab. ||). A factor 
c PT which can be used additionally to correct the measured cross section to the full phase 
space da/dQ 2 (p T > GeV 2 ) = c PT da/dQ 2 (p T > 25 GeV 2 ) is given in Tab. [| The error 
on c PT accounts for uncertainties in the structure functions. Results are given at fixed Q 2 
values, therefore the measured bin averaged cross section has been corrected by the factor 
Cbin = da d q2® V A AQ 2 ° r which is given in Tab. ^|. Differences in center of mass energy for 
different running periods are taken into account and the results are given for s = 90200 
GeV 2 . 

Agreement with the Standard Model prediction is observed. Our result for NC is also 
in agreement with the differential cross section measurements based on essentially the same 
data using the electron variables presented in |J with slightly different selection criteria. 
Above Q 2 = 5000 GeV 2 the cross sections for neutral and charged current processes are of 
comparable size. For the e + p data the charged current cross section stays below the neutral 
current cross section over the whole Q 2 range, due to the different y— dependence of the 
quark contributions to the cross section (cf. eq. (3,4), 0). 

The CC cross section errors are dominated by the statistical errors whereas for the NC 
samples the main contribution to the errors comes from the uncertainty in the hadronic 
energy scale. Other sources of systematic errors are uncertainties in the determination of 
the trigger and selection efficiencies as well as the luminosity measurement. 

For the neutral current cross section we show in Fig. 3 the Standard Model expectation 
for 7, Z° exchange and assuming photon exchange only. This illustrates that the present 
integrated luminosity does not yet allow the distinction between the pure photon exchange 
and the total electroweak cross section in NC reactions. 
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Figure 3: The unfolded differential cross sections da/dQ 2 are shown for charged and neutral current 
interactions with p? > 25 GeV for e~p collisions in a) and for e + p collisions in b). The inner error 
bars contain statistical, the full error bars include also the systematic errors added in quadrature. 
Boundaries of each bin are denoted by the vertical bars on the top of the plots. The solid lines 
indicate the Standard Model predictions for neutral current interactions and the dotted lines for 
charged current interactions. Dashed lines show the prediction for NC from photon exchange alone. 
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Table 1: Event statistics, differential, and total cross sections for px > 25 GeV are listed for NC 
and CC processes in e + p and e~p running. The first error is of statistical nature and obtained from 
the diagonal element of the covariance matrix, the second accounts for systematic uncertainties. 
The total cross sections are compared with the Standard Model prediction for m\y = 80.2 GeV 
and using the MRSA' parameterization for the proton parton densities. 



From the CC sample we obtain using the extended maximum likelihood fit a propagator 
mass of 

m w (e-p) = 78 t^ 1 ±j GeV , m w (e + p) = 97 ±{ 8 5 t 5 w GeV, 

where the first errors are of statistical and the second of systematic origin. The statistical 
uncertainties reflect the different sensitivities of the e + p and e~p cross sections to the prop- 
agator mass. The main source of systematic errors is the hadronic energy scale. The masses 
for both charge states are compatible and a combined fit results in 

m w = 84 ±1 + _\ GeV. 

Exploiting the Q 2 shape alone gives a consistent but less precise value for the propagator 
mass. 

The shape and cross section analysis exhibits that both charged current processes are 
consistent with the exchange of a If boson which has equal mass for both charge states. 
This confirms previous results from HERA on the W-m&ss 0, || with improved precision. 
In addition the mass of the intermediate W boson is in agreement with the mass determined 



from direct production of W bosons at pp-colliders p0 



In this analysis we use a larger angular coverage for the tracks entering the vertex re- 
construction compared to our previous analysis where the polar angular acceptance started 
at 15° ||. On the basis of this extended acceptance we update the total CC cross sections 
published in H. These cross sections as well as the ones for NC are collected in Tab. III. 
The Standard Model expectations are calculated as before using the HERACLES generator, 
where the error on the expectation originates in the choice of the proton structure function. 



We present also the ratios of total NC and CC cross sections R = unc/^cc- Due to the 
same treatment of neutral and charged current events many systematic errors cancel in this 
ratio. For p T > 25 GeV we obtain 



R+(e + p) = 12.8 ± 2.0 ± 1.0 



Br{e~ p) = 7.2 ± 1.6 ±0.7. 



11^ 



The result for R is in good agreement with the theoretical prediction of Bg M = 7. 2 ±0.1. We 
observe a two standard deviation discrepancy for R + from the prediction Bg M = 17.7 ±0.9. 

The residual systematic uncertainty is still caused by the hadronic energy scale due to the 
different cross section slope of NC and CC near the selection cut of px,h = 25 GeV. Selecting 
only events above Q 2 = 5000 GeV 2 , where the cross section slopes are nearly identical reduces 
the systematic uncertainty further. For such a cut we obtain R + = 2.0 ± 1.4 ± 0.07 and 
R~ = 1.4 ±0.8 ±0.05, while the expected values are R + = 4.17 and R~ = 1.34. With higher 
integrated luminosity these cross section ratios will be sensitive probes for precision tests of 
the Standard Model JHJ . 
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1.03 ±0.01 
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Table 2: For each Q 2 bin the bin boundaries, the applied correction factor for the finite bin size 
Chin and the correction factors to obtain a Born cross section &Born and to correct for the full phase 
space c PT are given. 



6 Conclusion 

The Q 2 dependence of the neutral and charged current processes have been measured for e~ 
and e + beams colliding with protons for transverse momenta of the scattered lepton above 
25 GeV. The analysis is based on the properties of the hadronic final state in both cases. The 
differential cross sections, as well as the total cross sections and their ratios are in agreement 
with the Standard Model expectations. At large momentum transfers Q 2 > 5000 GeV 2 



we observe that the weak and electroweak forces have the same strength. The propagator 
masses measured in charged current scattering are consistent with being the same for both 
charges. The average is my/ = 84 GeV and is consistent with the on-shell measurement 
of the W boson mass. 
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